function WT = WeightWT_1(q,par);
% WEIGHT_WT_1.M        Version I.K. - 22/12/06
% This function is called by the objective function OBJ_FCT_GMM_UNI_1
% Notations : Hamilton, ch 14

% Now at m = 0 out of q
gama0T = GamaM_1(0,par);
ST = gama0T;
m = 1;
while m <= q        % implementing eq 14.1.19 
    % Now at m out of q
%     m
    gamaMT = GamaM_1(m,par);
    ST = ST + (1-m/(q+1))*(gamaMT+gamaMT');
    m = m+1;
end
WT = inv(ST);
